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Abstract. The energy dependence of the local V and CP violation in Au+Au and Cu+Cu collisions in a large 
energy range is estimated within a simple phenomenological model. It is expected that at LHC the chiral magnetic 
effect will be about 20 times weaker than at RHIC. At lower energy range, covered by the low-energy scan at 
RHIC and future NICA/FAIR facilities, the created magnetic field strength and energy density of deconfined 
matter are rather high providing necessary conditions for the chiral magnetic effect. However, the particular 
model for the chiral magnetic effect predicts that this effect should vanish sharply at energy somewhere above 
the top SPS one. To elucidate CME background effects the Hadron-String-Dynamics (HSD) transport model 
including electromagnetic fields is put forward. Importance of new planning experiments at LHC and for the 
low-energy RHIC scan program is emphasized. 



1 Introduction 



2 Phenomenological estimates of the CME 



Until now there is no direct experimental evidence for topo- 
logical effects in QCD while the existence of nontrivial 
topological configurations is a fundamental property of the 
gauge theory. Transitions between different topological states 
occur with the local violation of the P and CP symmetry. 
Due to chiral anomaly the interplay of these topological 
configurations with (chiral) quarks results in asymmetry 
between left- and right-handed quarks. Such chiral asym- 
metry coupled to a strong magnetic field, created by col- 
liding nuclei perpendicularly to the reaction plane, induces 
a current of electric charge along the direction of a mag- 
netic field thereby separating particles of opposite charges 
with respect to the reaction plane. Thus, as was argued 
in Refs. B1I2I3I4I5I the topological effects in QCD may 
be observed in heavy ion collisions directly in the pres- 
ence of very intense external electromagnetic fields due 
to the "Chiral Magnetic Effect" (CME) as a manifesta- 
tion of spontaneous violation of the CP symmetry. First 
experimental evidence for the CME identified via the ob- 
served charge separation effect with respect to the reac- 
tion plane was presented by the STAR Collaboration at 
RHIC (6). In this paper we analyze the STAR data in a 
simple phenomenological way to estimate a possibility ob- 
serving the CME in the larger energy range, from the LHC 
to FAIR/NICA energies. We also make a step toward a 
dynamical estimate of the CME background based on the 
nonequalibrium Hadron-String-Dynamics (HSD) microscop- 
ical transport approach Q including electro-magnetic field. 



2.1 Model assumptions 

Qualitatively the CME may be estimated as follows. 

We consider the saturation momentum Q s as a charac- 
teristic scale of the process 0], so the transverse momen- 
tum of particles p, ~ Q s . Then the total transverse energy 
per unit rapidity at mid-rapidity deposited at the formation 
of hot matter is expressed through the overlapping surface 
of two colliding nuclei in the transverse plane S 



dy 



e-V = e-Az-S = Q s 



(Q]S) 



(i) 



Here the energy density and longitudinal size Az — At 
1 / Q s are taken in order of magnitude as follows e ~ Q\ 
and Az ~ 1/ Q s - The value in brackets is proportional to the 
hadron multiplicity (Q 2 S S) ~ dNhadrows/dy, thus we get 



dy 



dN hli . 



dy 



(2) 



which simply tells us that the total transverse energy per 
unit rapidity is equal to the total number of produced hadrons 
per unit rapidity times their average transverse momentum. 

The topological charge (winding number) generated in 
random motion during the time Tg when the magnetic field 
is present may be estimated as 



(3) 
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where / $ is the sphaleron transition rate. Sphalerons de- 
scribe a one dimensional random walk in the topological 
number space. At weak coupling, r$ ~ a 5 s T 4 . Little is 
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known about the rate of topological charge diffusion at 
strong coupling A = g 2 N c ; however in SUSY Yang-Mills 
theory the diffusion rate of topological charge can be eval- 
uated explicitly and is given by r$ ~ A 2 T 4 . 

The initial temperature To of the produced matter at 
time t — 1 IQ S is proportional to the saturation momen- 
tum Q s , Tq = c Q s . Neglecting the expansion time and the 
corresponding time dependence of the temperature we can 
roughly get 



= Ve 2 



^Vhadrons 



\ dy 



Vg7 



(4) 



Since sizable sphaleron transitions occur only in the de- 
confined phase, the time tb in Eq. (HJ is really the smallest 
lifetime between the strong magnetic field ?b one and the 
lifetime of deconfined matter T f : 



t b = min{r B ,r e j 



(5) 



The measured charged particle asymmetry is associ- 
ated with the averaged correlator a by the following rela- 
tion El: 

<cos(^ a + ^-2>P R p)>= (6) 
= (cos(i/r ff + ifr fi - 2i// c ~))/v 2 , c = v\, a v\,B - a a af3 , 

where Wrp is the azimuthal angle of the reaction plane de- 
fined by the beam axis and the line joining the centers of 
colliding nuclei and averaging in © is carried out over the 
whole event ensemble. The second equality in © corre- 
sponds to azimuthal measurements with respect to parti- 
cle of type c extracted from three-body correlation analy- 
sis J8|, v\ and vi are the directed and elliptic flow parame- 
ters, respectively. According to Ref. 0] an average corre- 
lator a = ^a a aB is related to the topological charge, n w , 
as 

a . (7) 

«Nhadrons/«I/ 

This estimate neglects absorption and rescattering in dense 
matter responsible for the difference of magnitudes between 
the same and opposite charge correlations and so is very 
rough. Nevertheless, for the purposes of our estimate, the 
use of Eqs. (|4|i and (0 yields the relation 



a 



h adrons /dy 
on which our consideration is based. 



(8) 



2.2 Energy dependence of the CME 

Let us now come to numerical estimates. Using the energy 
dependence Q 2 ~ s]lt ~ oWhadrons /dy deduced from RHIC 
and HERA |9 10 1 data we finally get for the energy depen- 
dence 




1/16 



(9) 



Thus, the direct energy dependence is comparatively weak. 
Results of dynamical heavy-ion calculations of the mag- 
netic field at the central point of the transverse overlap- 
ping region of colliding nuclei and energy density of cre- 
ated particles are presented in Figs. Q] and |2l respectively. 
Here for a field estimate we follow Ref. [ 1 1| basing on the 
UrQMD model B12I13II and applying the electromagnetic 
Lienard-Wiechert potential with the retardation condition 
for the magnetic field. As is seen, at the impact parameter 
b — 10 fm the maximal strength of the dominant magnetic 
field component B y (being perpendicular to the reaction 
plane) is decreased in Au+Au collisions by the factor of 
about 10, when one proceeds from y/s^N =200 GeV to 
Ei a b =11 GeV, while for the created particle energy den- 
sity e in the central box this factor is 250 i.e. noticeably 
higher. 



b=10fm 



— E, =10 A GeV 

lab 

- - E, =60 A GeV 

lab 

. - . E, =160 A GeV 

lab 




t [fm/c] 



Fig. 1. The time evolution of the magnetic field strength eB y at the 
central region in Au+Au collisions with the impact parameter b = 
10 fm for different bombarding energies. Calculations are carried 
out within the UrQMD model 1121131 (for a detail see (TTJ)- 



To use Eq. © we need to identify the impact param- 
eter, saturation momentum and multiplicity at a specific 
centrality. These can be found in Ref. [ 10] where the Glauber 
calculations were done. As a reference point we choose 
b =10 fm in our subsequent consideration. 

The measured value of {cos{i/i a + i//g - 2 'Frp)) for the 
same charge particles from Au+Au ( y^j^v = 

200 GeV) 

collisions at the impact parameter b =10 fm (40-50% cen- 
trality interval) is -(0.312 ± 0.027) • 10~ 3 0. Appropriate 
number for ysjviv = 62 GeV seems to be a little bit larger 
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Fig. 2. The time evolution of the energy density e of created 
particles in the Lorentz-contracted box with the 2 fm side at the 
central point of overlapping region. The impact parameter b = 10 
fm. 



but for Cu+Cu collisions the effect is definitely stronger [ 6 1 . 
Thus, ignoring any final state interactions with medium, 
assuming a a — og = a and neglecting the directed flow 
via = vib = we get from Eq. © a 2 exp — 0.31 • 10~ 3 for the 
maximal RHIC energy. Using numbers for the =200 
GeV reference case, from Eq. (0 we may quantify the CP 
violation effect by the correlator 



(10) 



The normalization constant Ka u can be tunned at the ref- 
erence energy y/lim =200 GeV from the inverse relation 
and experimental value a exp at this energy for b =10 fm 



Ka„ — 



exp 



(200) 



1/8 



r B (200) 



(11) 



The lifetime tb may be defined as the time during which 
the magnetic field is above the critical value needed to sup- 
port a fermion Landau level on the domain wall eS cn - ( = 
2n/Sd, where Sd is the domain wall area. Since the size 
of the domain wall is not reliably known, it is hard to pin 
down the number, but it should be of the order of m 2 . Hon- 
estly, we have to treat it as a free parameter. 

Indeed the size of the topological defect (say, a sphaleron) 
in the region between T c and 2T C is very uncertain. At 
weak coupling, the size is determined by the magnetic screen- 
ing mass and it is ~ l/(a s T). If one plugs a s «0.5 and T = 



Table 1. Estimated parameters for the CP violation effect in 
Au+Au collisions at centrality (40-50)% with the critical field 
eB crit = 0.2 ml. 



V^wv GeV 


1/16 


fa, fm/c 


t € , fm/c 


a 2 


4.5 ■ 10 3 


2.86 


0.018 


>1 


0.016- 10- 4 


200 


1.94 


0.24 


>2 


0.31-10 3 


130 


1.84 


0.33 


-2.3 


0.45-10- 3 


62 


1.68 


0.62 


-2.2 


0.93- 10~ 3 


17.9 


1.43 


1.41 


-2. 


2.48- 10~ 3 


11. 


1.35 


1.66 


- 1.9 


3.10-10- 3 


4.7 


1.21 


0. 


0. 


0. 



200 MeV, the size is of about 2 fm and then the threshold 
field is very small eB y ~ (a s T) 2 ~ 0.2 m 2 . 

On the other hand, we know that between T c and 2T C 
the magnetic screening mass which determines the size of 
the sphaleron is not small as expected from the pertur- 
bative theory, a s T, but from the lattice it is numerically 
large till about 5T C . This would increase the threshold to 
20 m 2 ,, however the relation between magnetic mass and 
the sphaleron size is valid only as long as the coupling is 
weak. 

All we can say it is perhaps in between (0.2-20) m\. 
Eventually lattice QCD calculations may clear this up. 

The upper bound on the magnetic strength eB cr ij — 
20 m 2 , results in = even for the RHIC energy and 
therefore in this case the CME should not be observable 
at all in this energy range. The time evolution of the mag- 
netic field and energy density, e, of newly created hadrons 
are presented in Figs. [TJ and [2] The extracted values of 
tb defined by the constraints eB y > 0.2 m 2 , and t £ (e > 
1 GeV/fm 3 ) are summed in Tabl. [TJ For the reference en- 
ergy and the minimal magnetic field constraint we have 



2.52 • 10" 



(12) 



If lifetimes are known for all energies one can estimate the 
CP violation effect through the a 2 excitation function. 

^From the first glimpse as follows from Tabl. [1] in the 
case of eB cr u = 0-2 m 2 , the interaction time tb is defined 
solely by evolution of the magnetic field since r# < t e 
whereas t e » 2 fm independent of yjs^N- The expected 
CME for Au+Au at b = 10 fm (see the last column in 
Tabl. [TJ monotonously increases when yjs^N goes down 
but then sharply vanishes exhibiting a shallow maximum 
in the range between near the top SPS and NICA energies. 
The position of CME maximum and its magnitude depend 
on the cut level which just defines fg- The decrease of the 
eBy bound to 0.02 m 2 , shifts the maximum toward lower 
energy -\[snn and enhances its magnitude. In an opposite 
limit when results are extrapolated to the LHC energy, the 
CME falls down by a factor of about 20 with respect to 
the RHIC energy. This result is quite understandable. The 
CME is mainly defined by the relaxation time of the mag- 
netic field which is concentrated in the Lorentz contracted 
nuclear region ~ 2R/y. Therefore, the CME inversely pro- 
portional to the colliding energy, - 1 / yfsNN and proceed- 
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ing from the RHIC to LHC energy we roughly get the sup- 
pression factor about 4.5/0.2~ 22. 

There is one worrying point here. Proceeding from 
yj/viv =200 to 62 GeV the predicted value of a 2 for b-\Q 
fm increases in three times though not more 20% growth 
has been observed in these collisions in the recent exper- 
iment [6]. This essential disagreement cannot be removed 
by a simple variation of eB cr j,. One may try to explain this 
correlator overestimation at yjs^N =62 GeV by an irrel- 
evant choice of the energy dependence of multiplicity in 
Eq. ( TTOb . For the correlator ratio at these two energies we 
have 

a 2 (200) _ T fl (200) / 62 \ 1/8 
a 2 (62) ~ r B (62) \200/ 

= 0.387(0.31/ ~ 0.72. (13) 

where we use lifetime values from Tabl.[T]and experimen- 
tal values for correlators |6|, B = 1/8. As follows from 
Eq. ( TOT l. to explain the experiment the exponent should be 
negative, B < 0. Therefore, the fast growth of t# with the 
energy decrease cannot be compensated by uncertainty in 
energy dependence of the correlator a. 

Uncertainty in the choice of the impact parameter does 
not help us to solve this issue as well. As is seen in Fig. [3] 
the maximal strength of the magnetic field slightly increases 
with the impact parameter but the lifetime of strong mag- 
netic field tb practically is independent of the impact pa- 
rameter (for example, look at the level of eB crll = 0.2 m 2 ). 
It will be used below as an estimate of the ^-dependence 
of the CME which is roughly proportional to b 0. 
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Fig. 3. The time evolution of the magnetic field strength eBy at 
the central region in Au-Au (200 GeV) collisions for different 
impact parameters. 

As follows from Eq. 0131 . to eliminate model discrep- 
ancy between two RHIC energies the lifetimes should be 
comparatively close to each other, Tg(62) w 1.2tb(200). 
One may try to fit this ratio by the variation of eB cra . As is 
seen from Tabl.[T]ateB cn 7 = 0.2 to 2 we have tb(62)/tb(200) ~ 
3.0. The needed time ratio, 1.2, is achieved at eB cr u ~ 
1.05 m 2 , very close to the maximal magnetic field of ~ 
1.2 m 2 for ys/wv=62 GeV (see FigQ]). It is evident that 
if the low-bound field eB cr n coincides with the maximal 
one the CME should not be observable. For more accu- 
rate tuning new experimental data at yfsim < 62 GeV are 



very needed. The decrease of the critical field till eB cr u ~ 
0.01 m 2 makes t# larger; however, the condition fg < r £ 
which could improved situation for tb is not yet achieved 
at ys'Aw =62 GeV. Unfortunately, the influence of this 
lower bound of the eB y insufficient for manifestation of 
this effect. 

The CME assumes that there are soft equilibrated quarks- 
gluons. But it is hard to think that there is an quark-gluon 
equilibration in the very early state and the concept of the 
initial time is introduced Q. It is used to suppose that an 
equilibrium state is achieved at the expansion stage which 
starts just after passing the maximum in the energy den- 
sity evolution. Thus, the t £ value should be count off from 
this time moment ?,( ysiviv) rather than from zero. This cor- 
rection makes shorter the lifetime by ?, i£ =0.08 and ~0.32 
fm for y7/vjv =200 and 62 GeV, respectively, as it follows 
from Fig. [2] Thus, the lifetime ratio at two energies dis- 
cussed is getting less, t b (62)/t b (200) * (0.62-0.32)/(0.24- 
0.08) ~ 2 instead of the mentioned «3.0. This softens the 
dependence of a 2 (200)/a 2 (62) on -\/Jnn but it still is not 
sufficient. 

According to Ref. [3| the initial time is estimated as 
ti ~ l/6.s, where Q s ~ as noted above. Normalizing 
the initial time at =200 GeV to t ue =0.08 fm we get 

Tb =0. 16 fm for 62 GeV. This shift is too short to influence 
essentially on the tb(20Q)/tb(62) ratio. With reference to 
the Lorentz contraction, other scenario maybe proposed as 
ti ~ 1/ Vsaw 131 141 . Under similar normalization we arrive 
at t, a:0.27 fm which is rather close to f; >£ at 62 GeV in the 
above estimate and therefore also is unsatisfactory. 

Now let us try to combine two last scenarios by chang- 
ing the low bound of the magnetic field eB cr u and simul- 
taneously shifting the appropriate interaction time to the 
initial moment defined by the energy density of created 
quarks, tb - TB(eB cr i,) - tj S . It is turned out that the equal- 
ity ( TT3l l is satisfied if eB cr u ~ 0.7 m 2 ,. Using the value of 
Tg(200) obtained in this analysis one can recalculate the 
coefficient in Eq. (fTOt . K - 6.05 • 10~ 3 , and therefore 
find the correlator a at any energy. In principle, similar 
analysis may be repeated for other impact parameters to 
consider the ^-dependence of the CME. As was shown in 
Refs. 01141 the CME roughly is linear in b/R. Taking this 
as a hypothesis we evaluate the centrality dependence of 
the CME fitting this line to points b - 10 fm (or centrality 
(40 - 50)%) to be estimated in our model and b - where 
the CME is zero. The results are presented in Fig. [4] for 
Au+Au collisions at three energies. 

As it is seen the calculated lines quite reasonably re- 
produce the measured points of azimuthal asymmetry of 
charge separation for Au+Au collisions at yjs^N =200 and 
62 GeV. The chosen value of eB cr j, = 0.7 m 2 results in ab- 
sence of the CME at the top SPS energy because the crit- 
ical magnetic field practically coincides with the maximal 
field at this bombarding energy (see Fig. [TJ. The CME at 
the LHC energy is expected to be less by a factor of about 
20 as compared to that at the RHIC energy. Note that at the 
LHC energy we applied a simplified semi-analytical model 
ifTTl for magnetic field creation and assumed f, £ = 0. Thus, 
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Fig. 4. Centrality dependence of the CME. Experimental points 
for Au+Au and Cu+Cu collisions are from |6|. The dotted line is 
our prediction for Au+Au collisions at the LHC energy. 



we consider this LHC estimate as an upper limit for the 
CME. 



2.3 System-size dependence of the CME 

Similar analysis may be repeated for Cu+Cu collisions bas- 
ing on available RHIC measurements at two collision en- 
ergies. Here one remark is in order. An enhancement of the 
CME in Cu+Cu collisions with respect to Au+Au ones was 
seen experimentally at the same centrality |6| but not at the 
same impact parameter. As follows from the Glauber cal- 
culations [ 15], the impact parameter b=10 fm for gold re- 
actions corresponds to centrality (40-50)% while the same 
centrality for copper collisions matches b=4.2 fm. As is 
seen in Fig(5]the time distributions of the magnetic field 
and energy density for Cu+Cu collisions look very simi- 
lar to that for Au+Au ones but lifetimes, both ?b and r £ , 
are shorter in the Cu+Cu case. The extracted lifetimes and 
other characteristics for <?fi cnr = 0.2m 2 (Kc u = 6.34 • 10~ 3 ) 
are presented in Tabl. [2] We meet again the same prob- 
lem: one should compensate a too strong energy depen- 
dence of the model correlators by the proper definition of 
lifetimes. Defining the lifetime in the same manner as for 
Au+Au collisions tb - fB(eB c ,n) - t ue the lifetime ratio 
at eB cr i t = 0.3 m 2 for Cu+Cu collisions tb(62)/tb(200) = 
(0.35 - 0.21)/(0.14 - 0.06) is turned out to be very close to 
experimental one. In this case K Cu = 11.9 - 10 3 . In the lin- 
ear approximation with the reference point at b =4.2 fm, 
one may draw the centrality dependence of the CME for 
Cu+Cu collisions shown also in Fig|4]which is in a reason- 
able agreement with the experiment. Note that eB crll =0.3 
m% which is slightly above the maximal magnetic field at 
y/s NN =62 GeV implies that the CME for Cu+Cu colli- 
sions will not be observable even at the top SPS energy. 

From dimensionality arguments the system-size depen- 
dence of the chiral magnetic effect (at the same all other 
conditions) would be expected to be defined by the surface 
S = S a{£>) of an "almond"-like transverse area of over- 
lapping nuclei since both the high magnetic field and de- 
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Fig. 5. The time evolution of the magnetic field strength eB,j and 
deconfined energy density e at the central region in Cu+Cu col- 
lisions at a/sjw =200 and 62 GeV for centrality (40-50)% (b = 
4.2 fm). 



Table 2. Estimated parameters for the CP violation effect for 
Cu+Cu (b=4.2 fm) collisions with eB cril = 0.2 m\. 



VW GeV 


1/16 
S NN 


r fl , fm/c 


r £ , fm/c 


a 2 


200 


1.94 


0.15 


>2. 


0.49- 10~ 3 


62 


1.68 


0.37 


>2. 


1.40- 10" 3 



confined matter are needed for this effect. The magnetic 
field was evaluated in the the center of the overlapping re- 
gion but as was shown in Ref. ifTTIl the studied eB y com- 
ponent is quite homogeneous along y of this "almond". 
Above, say in the -^sJv/v"-analysis of Au+Au collisions, the 
area ,5>Au(b=10 fm) does not depend on the energy and it 
is involved into the coefficient K, Eq. ([Til l. Using for "al- 
mond" area a rough estimate as for sharp boundary over- 
lapping discs of radius R - roA 1 ^ 3 , namely S = S aQ>) = 
n^R 2 -(b/2) 2 (R - b/2), we have S C u(b = 4.2) /S Au (b = 
10) ~ 1.65 which seems to be consistent with experimen- 
tal ratio of the CME at these two points. Thus, if one ne- 
glects the lifetime on the system-size and keep only the 
renormalized coefficient Kq u = K^u Scu(b = 4.2)/5a u (^ 
1 j 1.65,^11 f° r the Au curve at -\/snn =200 GeV in 
Fig|4]we are able to reproduce rather well the results for 
Cu+Cu (200 GeV). Note that this Au curve was obtained 
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for eB cr i t = 0.7 and non-zero initial time f, £ . Thus, con- 
ditions of the compared events are not identical and this 
success cannot be repeated for Cu+Cu (62 GeV) collisions 
where maximum of the magnetic field is lower than 0.7 m^. 
The above rough estimate essentially depends on the ro pa- 
rameter. From our analysis of the centrality dependence of 
the CME at =200 GeV at eB cr u = 0.2 m\ we have 

K Cu (b = 4.2)/K Au {b = 10) = 6.34 • 10- 3 /2.52 . 10" 3 * 2.5 
which strongly differs from the above estimate. Therefore, 
the Cu enhancement effect is not only a geometric one. In 
this case if the reference point in centrality is considered, 
the energy dependence should be taking into account. One 
should note that in a general case the saturation momen- 
tum Q% is weakly A-dependent as well what was neglected 
in our consideration. 



3 Towards a kinetic approach to the CME 
background 

The discussed CME signal, the electric charge asymme- 
try with respect to the reaction plane, may originate not 
only from the spontaneous local CP violation but also be 
simulated by other possible effects. In this respect it is im- 
portant to consider the CME background. We shall do that 
considering a full evolution of nucleus-nucleus collisions 
in terms of the HSD transport model [7 1 but including for- 
mation of electromagnetic field as well as its evolution and 
impact on particle propagation. 

Generalized on-shell transport equations for strongly 
interacting particles in the presence of electromagnetic fields 
can be written as . 



ij t + ( v p u ) v r - (V, U + qy x (V x A)) V p } /(r, p, t) 



- Icoll(f,fl, — /iv) 



(14) 



which are supplemented by equations for magnetic and 
electric fields 

<9A 

B = VxA, E = -V<5- — . (15) 

ot 

The hadronic mean-field U ~ Re(£ re ')/2po here has an 
additional term from the 4-vector electromagnetic poten- 
tial A = {0, A). The general solution of the wave equa- 
tions (fT~5b with charge distribution p(r, f) = qn and current 
j(r, f) = q\ are 

1 r P (r',t')6(t-t'-\r-r'\/c) , ,. lf 



4tt J |r-r'| 

for the electromagnetic potential @(r, t) and 

> , A 1 f j(r',Q5(f-f'-|r-r'|/ c ) j3 
A(r,r)=-J j— drdt(\l) 

for the vector potential. For moving point-like charges, one 
gets 

p(r, i) = q 6(r - r(f)); j(r, i) = q v(f) 6(r - r(t)). (18) 



In this case Eq. (fTFt leads to the above used retarded Lienard- 
Wiechert potentials generated by every moving charged 
particle and acting at the point R = r - r' at the moment 

/CD. 
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Fig. 6. Distribution of the magnetic field strength eB,j in the y = 
plane at / =0.05 (upper panel), 0.2 (in the middle) and 0.55 (lower 
panel) fm/c for Au+Au collisions at ^]s NN =200 and b =10.2 
fm. The collision geometry is projected on x - z plane by points 
corresponding to a particular spectator position. Curves (and their 
projections) are levels of the constant eB r 



One should note that the off-shell HSD transport ap- 
proach is based not on the Boltzmann-like transport equa- 
tion (TBI) but rather on the off-shell Kadanoff-Baym equa- 
tions [ 16 1 having similar general structure. The set of equa- 
tions was solved in a quasiparticle approximation by us- 
ing the Monte-Carlo parallel ensemble method. To find the 
magnetic field a space grid was used. In a lattice point of 
this grid the retarded vector potential is calculated by nu- 
merical differentiation. The field inside a cell is approxi- 
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mated by that at the nearest grid point. To avoid singulari- 
ties and self-interaction effects, particles within a given cell 
are excluded from procedure of the field calculation. 



AuAu, Vs AW = 20OGeV, b=10.2fm, t=0.05fm/c 
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Fig. 7. The AC — z projection of the energy density e and magnetic 
field strength eB tJ (contour lines) in the y = plane for three cases 
presented in Fig[6] 



An evolution snapshot of the magnetic field B y {x, y = 
0, z, t) (in units of m 2 n ) formed in Au+Au (200 GeV) pe- 
ripheral (b =10.2 fm) collisions are given in Fig|6]for three 
time moments t =0.05, 0.20 and 0.55 fm/c. The collisional 
geometry is presented by a set of points every of which 
corresponds to a spectator nucleons. The field is not homo- 
geneous exhibiting a wide maximum over the transverse 
size of overlapping (participant) matter. Opposite rotation 
of the magnetic field along direction of two colliding nu- 
clei results in corresponding two minima from outer sides 
of spectator matter remnants. At expansion these remnants 
are moving away from each other. It is of interest that for 



the fixed x the width of Z distributions of the magnetic field 
does not changed much. This type of scaling has been as- 
sumed in the phenomenological model of the CME consid- 
ered above. The maximum in the magnetic field strongly 
correlates with that in the energy density of created par- 
ticles as is clearly seen from Figj7] Large local values of 
By and s reached in these Au+Au collisions provide nec- 
essary conditions for observation of signals of a possible 
parity violation. 



AuAu, VSj/u = 200GeV, b=10.2fm 




Fig. 8. The time evolution of the magnetic field eBy at the point 
where a pion is situated (shown by points). The curve corresponds 
to the field averaged over these points at every time moment. 



A particular characteristic is shown in Fig|8] The mag- 
netic field in the point where a real meson is situated is 
plotted as a function of time for every meson. The curve 
corresponds to the averaged field < B y /m^ > at the given 
time moment. It is seen that averaged magnetic field sharply 
goes down till B y < m 2 and then falls down quite slowly. 
The moment, when the expansion regime is changed, is of 
order of 2R/y cm ~ 0.2 fm/c. 



4 Discussion and conclusions 

Summarizing Sectf2] one should note that for heavy-ion 
collisions at y%5v > 1 1 GeV the magnetic field and en- 
ergy density of deconfined matter reach very high values 
which seem to be high enough for manifestation of the 
Chiral Magnetic Effect. However, these are only necessary 
conditions. To estimate a possible CME a particular model 
is needed. For the average correlator our qualitative predic- 
tion a 2 ~ s~^ 8 has a rather small exponent but nevertheless 
it is too strong to describe observable energy behavior of 
the CME. This model energy dependence can be reconcile 
with experiment [6] by a detailed treatment of the lifetime 
taking into account both the time of being in a strong mag- 
netic field and time evolution of the energy density in the 
QGP phase. For the chosen parameters we are able to de- 
scribe data for Au+Au collisions on charge separation at 
two available energies. We predict that the effect will be 
much smaller at the LHC energy and will sharply disap- 
pear near the top energy of SPS. Experiments at coming 
into operation Large Hadron Collider and that planned at 
RHIC by the Beam Energy Scan program ifTTl are of great 
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interest since they will allow one to test the CME scenario 
and to infer the critical magnetic field eB cr i t governing by 
the spontaneous local CP violation. 

All the results presented above are obtained under as- 
sumption that evolution of the electromagnetic field scales 
as that in the central point of the overlapping region. The 
fall of By in time is caused be the fact that spectators are 
flying away. This assumption is getting worse for lower 
SPS energies what should be tested in more elaborated 
models. 

The experimentally observed CME enhancement for 
Cu+Cu collisions is related with the selection of different 
impact parameters for the same centrality at variation of 
a colliding system. However, it is not reduced to a purely 
geometrical effect. 

The problem of parity violation in strong interactions 
and the related CME are actively debated now. In particu- 
larly, the important finding has been reported recently [ 1 8 1 . 
To distinguish between the CME and effects due to certain 
two-particle correlations it was proposed to analyze the 
strength of charge separation effect and its azimuthal corre- 
lation with the reaction plane. It turned out that the observ- 
able effect for like-sign pions is the sum of a large negative 
in-plane correlator and small positive out-of-plane one, i.e. 
the correlations for same charge pairs are mainly in-plane. 
For unlike pions both components are positive and approx- 
imately equal to each other. All that is puzzling since the 
CME is expected only for out-of-plane events. A full dy- 
namical simulation is needed for heavy-ion collisions with 
inclusion of the magnetic field formation, its evolution and 
its impact on particle trajectories. First step in this direction 
has been made in Sec|3] Such approach allows one to cal- 
culate the CME background. It is of great interest that the 
discussed CME signal, the electric charge asymmetry with 
respect to the reaction plane, may originate not only from 
the spontaneous local CP violation but also be simulated by 
other possible effects. The search for alternative explana- 
tions and additional manifestations of local parity violation 
is underway 11191201211221231241251261271281291301311321331341 ). 
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It is important that the developed kinetic approach in prin- 
ciple allows one to simulate the Chiral Magnetic effect it- 
self. This work is in progress. 
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